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We use the numerically exact superposition T-matrix method to compute the optical cross sections and the Stokes 
scattering matrix for polydisperse mineral aerosols (modeled as homogeneous spheres) covered with a large number 
of much smaller soot particles. These results are compared with the Lorenz-Mie results for a uniform external 
mixture of mineral and soot aerosols. We show that the effect of soot particles adhering to large mineral particles 
can be to change the extinction and scattering cross sections and the asymmetry parameter quite substantially. The 
effect on the phase function and degree of linear polarization can be equally significant. © 2012 Optical Society of 
America 

OCIS codes: 290.0290, 290.1090, 290.5825, 290.5850. 


Mineral and soot particles are ubiquitous constituents of 
natural and artificial environments [1-5]. Both particle 
types can have quite complex morphologies [2—4] , there- 
by making the calculation of their scattering and absorp- 
tion properties highly nontrivial and necessitating the use 
of a combination of computer solvers of the macroscopic 
Maxwell equations with supplementary ranges of perfor- 
mance characteristics [6,7]. As if to aggravate the pro- 
blem even further, the morphologically complex mineral 
and soot particles can stick to each other and form var- 
ious types of compound aerosols. Whether the optical 
properties of a mixture of mineral and soot particles 
change upon aggregation depends on the strength of 
the electromagnetic interaction between the mineral and 
soot components, which can now be in direct physical 
contact instead of being widely separated. This interac- 
tion and the resulting optical effects are still insufficiently 
understood and poorly quantified, and so they constitute 
the main subject of this Letter. 

Depending on the specific type of environment, one 
can expect a varying range of morphologies of the indi- 
vidual mineral and soot particles. However, the main 
focus of this Letter is not on the optical effects of 
individual-particle morphology (see [6,8-12]) but rather 
on the potential effects of small soot particles adhering 
to the surfaces of larger mineral hosts. It is therefore ap- 
propriate to constrain the scope of this study by choosing 
the simplest morphology of the individual mineral and 
soot aerosols in the form of homogeneous spheres. This 
choice also allows us to use the highly efficient and nu- 
merically exact superposition T-matrix method (STMM) 
[13,14], 

Our model of a compound mineral-soot particle is illu- 
strated by the inset in Fig. 1_, wherein the surface of a 
relatively large mineral particle is covered by a number 
N s of small soot particles. This model differs from those 
analyzed in [10] in that the soot particles are assumed to 
be compact, thereby corresponding to the stage in their 
aging when fluffy soot fractals have already collapsed 
into densely packed aggregates with fractal dimensions 


approaching 3 [4,15]. The ratio of the radius of the miner- 
al host to that of the small soot particles is kept constant 
at R/r =10. Soot particles are placed on the surface of a 
large mineral host using a random-number generator and 
making sure that the particles do not overlap. The scat- 
tering and absorption characteristics of the compound 
aerosols are then averaged over the uniform orientation 
distribution of the resulting cluster as well as over a stan- 
dard power law distribution n{R) — constant x R~ 3 [8] 
with an effective radius R ( , n — 1 pm, effective variance 
u e£f = 0.05, and host radii in the range 0.661 pm 
<R< 1.439 pm. The resulting effective radius of the soot 
particles (r eff = 0.1 pm) is typical of the postcollapse 
stage [4,15]. The refractive index of the mineral hosts 
is 1.55 + i0.0003, while that of the soot particles is 
1.75 + 70.435. The wavelength is fixed at 0.62832 pm. The 
averaging over cluster orientations is performed using 
the efficient semianalytical procedure afforded by the 
STMM [13]. The numerical size averaging is based on 
the Gauss-Legendre quadrature formula with 100 nodes. 
The size averaging has a dual purpose of (i) better repre- 
senting polydisperse particle ensembles encountered in 
various applications and (ii) smoothing out artifacts in 
the form of sharp resonance features typical of monodis- 
perse scatterers [8]. 

The solid curves in Fig. 1 visualize the elements of the 
normalized Stokes scattering matrix relating the Stokes 
parameters of the incident quasi-monochromatic light 
and those of the scattered light: 
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where 0 e [0°, 180°] is the angle between the incidence 
and scattering directions; both sets of the Stokes para- 
meters are defined with respect to the scattering plane 


0146-9592/12/040704-03$15.00/0 


© 2012 Optical Society of America 


February 15, 2012 / Vol. 37, No. 4 / OPTICS LETTERS 705 



Fig. 1. Solid curves depict the elements of the size- and orientation-averaged Stokes scattering matrix for a large mineral particle 
covered randomly by N s small soot particles. The inset illustrates the morphology of the compound aerosols and shows a mineral 
host covered with N s = 49 soot particles. Dashed curves depict the Lorenz-Mie results for a composition-equivalent external mix- 
ture of mineral and soot aerosols. 


[8]. The (1,1) element of the scattering matrix is the con- 
ventional phase function normalized according to 

^ [ d0 sin 0Gq(0) = 1. (2) 

^ Jo 

The corresponding integral radiometric characteristics 
are listed in Table 1. These STMM results are paralleled 
by the results of Lorenz-Mie computations based on the 
assumption that each mineral host and the respective N s 
soot particles are all widely separated and scatter light 
independently of each other. 

The comparison of these two sets of results can be 
used to identify and quantify specific optical effects of 
aggregation. It is obvious that depending on A ? s . , these ef- 
fects can be quite strong. At least some of them can be 
interpreted qualitatively. 

First, increasing N s always increases the extinction 
and absorption cross sections, but this increase is slower 
in the case of the compound aerosols because the large 
mineral hosts shade almost half of the soot particles and 
thus preclude them from interacting with the incident 
field. Second, increasing the number of strongly absorb- 
ing soot particles always results in a decrease of the sin- 
gle-scattering albedo. It is interesting, however, that this 
effect is virtually the same for compound aerosols and 
the uniform external mixture of mineral and soot aero- 
sols. Third, the effects of increasing N s on the scattering 
cross sections and asymmetry parameters of the com- 


pound aerosols and of the uniform external mixture 
are exactly opposite. The asymmetry parameter of the 
external mixture decreases with A\ s . because the asymme- 
try parameter of the Rayleigh-like soot particles is much 
smaller than that of the large mineral hosts. In the case of 
compound aerosols, increasing N s causes the formation 
of a strongly absorbing coating around the weakly ab- 
sorbing mineral core, especially when soot particles 
cover a significant fraction of the host’s surface. The re- 
sulting phase function is increasingly controlled by the 
forward-scattering diffraction peak and the first external 
reflection, thereby becoming progressively anisotropic. 
This coating also suppresses the contribution of the host 
to the scattering cross section. Overall, Table 1 demon- 
strates that aggregation can result in a decrease of the 
optical cross sections by almost a factor of 2 as well as 
in significantly larger asymmetry parameters approach- 
ing 0.9. 

The left-hand panel of Fig. J_ demonstrates potentially 
very strong effects of aggregation on the phase function. 
For example, the solid and dashed blue curves exhibit a 
factor of 4.5 differences at side scattering angles and 1 
order of magnitude differences at backscattering angles. 
The middle panel of Fig. 1 is also indicative of potentially 
strong effects of aggregation on the degree of linear 
polarization for unpolarized incident light given by 
P(Q) — -b\ (0) /«] (0). In particular, the solid and 
dashed blue curves exhibit at backscattering angles a 
strong positive maximum and a deep negative minimum, 


Table 1. Ensemble-Averaged Extinction, Scattering, and Absorption Cross Sections (in //m 2 ) per Compound 
Aerosol Particle, Single-Scattering Albedo co, and Asymmetry Parameter g as Functions of the Number of Small 

Soot Particles Covering the Surface of a Large Mineral Host" 
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“The primed characteristics correspond to the case of a uniform external mixture of widely separated, independently scattering mineral and soot 
aerosols 
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respectively, which almost mirror each other with re- 
spect to the horizontal P(&) — 0 line. 

The ratio a 2 (0)/cfi(0) is identically equal to unity for 
the external mixture of spherical mineral and soot parti- 
cles but deviates from unity noticeably for compound 
aerosols, thereby signifying the nonsphericity of the lat- 
ter (see the right-hand upper panel of Fig. 1). The change 
of this ratio with N s is highly nonmonotonic. By the time 
N s reaches the value 249, soot particles form a quasi- 
uniform coating around the mineral host, thereby making 
the compound particle quasi-spherical. This appears to 
explain the smaller deviation of the blue curve from 
100% than those of the green and pink curves. Overall, 
the deviation of the ratio a 2 (0)/ai (0) from unity appears 
to be too weak to make this ratio useful in optical particle 
characterization. 

The ratio b 2 (&)/a l (&) appears to be even less sensi- 
tive to the presence of soot particles. Indeed, the right- 
hand bottom panel of Fig. 1 as well as the respective 
Lorenz-Mie results for the external mixture of mineral 
and soot aerosols (not shown) show rather weak 
changes with increasing N s . 

In summary, our numerically exact STMM and Lorenz- 
Mie results demonstrate that the effect of small soot par- 
ticles adhering to large mineral particles can be to change 
the optical cross sections and the asymmetry parameter 
quite substantially. The effect on the phase function and 
degree of linear polarization can be equally dramatic. 
There is no reason to expect that these conclusions will 
change significantly if one considers nonspherical rather 
than spherical mineral hosts, especially in the case of the 
integral radiometric characteristics. 

These conclusions are fundamentally different from 
those drawn in [16] on the basis of similar computations 
but performed for large mineral hosts covered with much 
smaller mineral grains of the same chemical composi- 
tion. There is no doubt that this difference is primarily 
caused by soot particles being much more absorbing 
than size-equivalent mineral grains. 

We expect that our results will be important in analyses 
of laboratory and remote-sensing data as well as in assess- 
ments of radiative aerosol effects on climate (e.g., 
[17-19]). The quantitative significance of our findings 
may vary depending on the type of environment and 
the specific application. The strong dependence of the ef- 
fects of aggregation on N s obviously indicates the urgent 
need for accurate measurements of the state of aggrega- 
tion of mineral and soot aerosols as well as of their relative 
number densities in various types of environment. 
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